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ABSTRACT 

We show that the epoch(s) of reionization when the average ionization fraction of the universe is 
about half can be determined by correlating Cosmic Microwave Background (CMB) temperature maps 
with 21-cm line maps at degree scales (/ ~ 100). During reionization peculiar motion of free electrons 
induces the Doppler anisotropy of the CMB, while density fluctuations of neutral hydrogen induce the 
21-cm line anisotropy. In our simplified model of inhomogeneous reionization, a positive correlation 
arises as the universe reionizes whereas a negative correlation arises as the universe recombines; thus, 
the sign of the correlation provides information on the reionization history which cannot be obtained by 
present means. The signal comes mainly from large scales (fc ~ 10“^ Mpc~^) where linear perturbation 
theory is still valid and complexity due to patchy reionization is averaged out. Since the Doppler signal 
comes from ionized regions and the 21 -cm comes from neutral ones, the correlation has a well defined 
peak(s) in redshift when the average ionization fraction of the universe is about half. Furthermore, 
the cross-correlation is much less sensitive to systematic errors, especially foreground emission, than 
the auto-correlation of 21 -cm lines: this is analogous to the temperature-polarization correlation of 
the CMB being more immune to systematic errors than the polarization-polarization. Therefore, we 
argue that the Doppler-21cm correlation provides a robust measurement of the 21-cm anisotropy, 
which can also be used as a diagnostic tool for detected signals in the 21 -cm data — detection of 
the cross-correlation provides the strongest confirmation that the detected signal is of cosmological 
origin. We show that the Square Kilometer Array can easily measure the predicted correlation signal 
for 1 year of survey observation. 

Subject headings: cosmic microwave background - cosmology: theory - diffuse radiation - galaxies: 
formation - intergalactic medium 


1. INTRODUCTION 

When and how was the universe reionized? This question is deeply connected to the physics of formation and 
evolution of the first generations of ionizing sources (stars or quasars or both) and the physical conditions in the 
interstellar and the intergalactic media in a high redshift universe. This field has been developed mostly theoretically 
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Each one of these methods probes different epochs and aspects of cosmic reionization: the Gunn-Peterson test is 
sensitive to a very small amount of residual neutral hydrogen present at the late stages of reionization ( 0 ^ 6 ), Lya- 
emitting galaxies and the wavelengths of the near infrared background probe the intermediate stages of reionization 
(7 ^ 15), the 21 -cm background probes the earlier stages where the majority of the intergalactic medium is still 

neutral (10 < z < 30), and the CMB polarization measures the column density of free electrons integrated over a 
broader redshift range (z ^ 20, say). Since different datasets are complementary, one expects that cross-correlations 
between them add more information than can be obtained by each dataset alone. For example, the information content 
in the CMB and the 21-cm background cannot be exploited fully until the cross-correlation is studied: if we just extract 
the power spectrum from each dataset, we do not exhaust the information content in the whole dataset because we 
are ignoring the cross-correlation between the two. The cross-correlation always reveals more information than can be 
obtained from the datasets individually unless the two are perfectly correlated (and Gaussian) or totally uncorrelated. 

Motivated by these considerations, we study the cross-correlation between the CMB temperature anisotropy and 
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the 21-cm background on large scales. We show that the CMB anisotropy from the Doppler effect and the 21-cm 
line background can be antz-correlated or correlated at degree scales {I ^ 100), and both the amplitude and the sign 
of the correlation tell us how rapidly the universe reionized or recombined, and locations of the correlation (or the 
anti-correlation) peak(s) in redshift space tell us when reionization or recombination happened. This information is 
difficult to extract from either the CMB or the 21-cm data alone. Our work is different from recent work on a similar 
subject bv iSalvaterra et ^ (|2nnfi(i . While they studied a similar cross-correlation on very small scales (~ arc-minutes), 
we focus on much larger scal es (^ degrees) where matter fluctuations are still linear and complexity due to patchy 
reionization is averaged out. iCooravl 1200411 studied higher-order correlations such as the bispectrum on arc-minute 
scales. For our case, however, fluctuations are expected to follow nearly Gaussian statistics on large scales, and thus 
one cannot obtain more information from higher-order statistics. He also studied the cross-correlation power spectrum 
of the CMB and projected 21-cm maps, and concluded that the signal would be too small to be detectable owing to 
the line-of-sight cancellation of the Doppler signal in the CMB. However, we show that cancellation can be partially 
avoided by cross-correlating the CMB map with 21-cm maps at different redshifts (tomography). Prospects for the 
Square Kilometer Array (SKA) to measure the cross-correlation signal on degree scales are shown to be promising. 

Throughout the paper, we use c = 1 and the following convention for the Fourier transformation: 

/(n,,7)=|^/ke-**^-(’'»-''), (1) 


where n is the directional cosine along the line of sight pointing toward the celestial sphere, ry is the conformal time, 
ri{z) = /p dt'/a{t') = dz'/H{z'), and ryo is the conformal time at present. Note that 


m - V{z) 


r dz' 

Jo H{z'Y 


( 2 ) 


which equals the comoving distance, r{z) = rjo — r](z), in flat geometry (with c = 1). Also, using Rayleigh’s formula 
one obtains 


/(n,ry) = 47r^(-i)' / —^/kjj[fc(ryo - T])]Yim{n)YiY{k). 

Im d \ 


( 3 ) 


The cosmological parameters are fixed at = 0.3, D;, = 0.046, Ha = 0.7, h = 0.7, and as = 0.85, and we assume a 
scale invariant initial power spectrum for matter perturbations. 

This paper is organized as follows. In §|2and|21we derive the analytic formula for the Doppler-21-cm correlation 
power spectrum. Equations US) or are the main result. We then present a physical picture of the correlation 
and describe properties of the correlation in detail. We also discuss the validity of our assumptions and possible 
effects of more realistic reionization scenarios. In §2|we discuss detectability of the correlation signal with SKA before 
concluding in § [S] 


2. 21-CM FLUCTUATIONS AND CMB DOPPLER ANISOTROPY 
2.1. 21-cm Signal 

Following the notation of IZaldarriaga. Furlanetto. fc Hernauisd ll2004H . we write the observed differential brightness 
temperature of the 21-cm emission line at A = 21 cm(l -|- z) in the direction of n as 


721(11,2) 


rr]o 

To{z) / dT]'W[r]{z) 

Jo 
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where W{ri{z) — rj') is a normalized (/.^ dxW{x) = 1) spectral response function of an instrument which is centered 
at ri{z) — ry' = 0, To{z) is a normalization factor given by 


To ( 2 ) ~ 23 mK 





and 

V^2i(n,7y) = a;ff(n,7y)[l-|-(5b(n, ry)] 1 
where Si, is the baryon density contrast. 
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is the ionized fraction contrast, Xh is the neutral fraction, and Xe = 1 — xr is the ionized fraction. Here, we have 
assumed that the spin temperature of neutral hy drogen, Ts. is much larg er than the CMB temperature, Tcmb- This 
assumption is valid soon after reionization begins llCiardi fc MadaiJl2003ll . 

To simplify the calculation, we assume that the spectral resolution of the instrument is much smaller than the 
features of the target signal i n redshift space. This is always a v e ry go od approximation. (For the effect of a 
relatively large bandwidth, see tZaldarriaga. Furlanetto. fc HernauistI ((20^.) Therefore, we set W{x) = 5^{x) to 
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obtain T 2 i(h, z) = To(z)'!/' 2 i[n, 77(z)]. To leading order in 5x and Sb, the spherical harmonic transform of T 2 i(n,z) is 
given by 

/ (Pk 

j^^h{z)(\ + f^x^)5b\, - Xe{z)5xv]a^^{k,z)Yi^{\i), ( 8 ) 

where a^^(k, z) is a transfer function for the 21-cm line, 

= To{z)D{z)ji[k{r]Q - r])], (9) 

D{z) is the growth factor of linear perturbations, /i = k-n, and / = dlnZ?/dlna. The factor (1 + /^^) takes account of 
the enhancement of the fluctuation amplitude due to the redshift-space distortion, the so-called “Kaiser effect” (Kaiser 
1987; see also Bharadwaj & Ali 2004 and Barkana & Loeb 2005). 

2.2. Doppler Signal 

The CMB temperature anisotropy from the Doppler effect is given by 

rvo rvo t d^k 

TD{n)=-TcmhJ d?7Te“”ri • Vf,(n,77) =-Tcmby d7?Te“”n • y (10) 

where Tcmb = 2.725 K is the present-day CMB temperature, x = dx/dp, t{p) = ot Jq dp'ne{p') is the Thomson 
scattering optical depth, and Vbk is the peculiar velocity of baryons. In deriving the above formula, we have neglected 
the fluctuation of ionized fraction, and electron number density, Sein,p), since their contributions to the 

cross-correlation would be higher order corrections (the effect due to is called the Ostriker-Vishniac effect; e.g., 
Ostriker & Vishniac 1986), and such a correction is negligible for linear fluctuations on the large scales we consider here. 
Note that the negative sign ensures that we see a blueshift, Tb){h) > 0, when baryons are moving toward us, n-<0. 
The peculiar velocity is related to the density contrast via the continuity equation for baryons, Vf,k = —i(}i/k‘^)5b\^D. 
One obtains 

r.(ft) = <»-.>. (11, 

The spherical harmonic transform of T]j{n, z) is then given by 

/ d^k 

^6b^aPik)Y;^{k), (12) 

where af‘{k) is a transfer function for the Doppler effect, 

aP(fc) = ^^y dpDTe~'^-^ji[k{po-p)]. (13) 


3. DOPPLER-21-CM CORRELATION 
3.1. Generic Formula 

Given the spherical harmonic coefficients just derived for the 21-cm line (Eq. |S|) and the Doppler anisotropy 
(Eq. [T^ b one can calculate the cross-correlation power spectrum, Cf^ ^, exactly as 


2 

= = - k^dk[xHiz){l +f{p.^))Pssik)-Xe{z)Pxs{k)]a‘f^{k,z)a[’{k) 
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where we have defined the matter power spectrum, Pss{k), as ((^kf^k') — (27r)^(5(k—k')P 55 (A;), and the cross-correlation 
power spectrum between ionized fraction and density Pxs{k), as (i5a;k(5((,) = (27r)^(5(k — k')Pxs{k). In the last line of 
equation itTH) . we have used f{p?‘) = 1/3 for a matter-dominated universe. Note that 5 used in these power spectra is 
the density contrast of total matter, 6, as baryons trace total matter perturbations, 6b = 6, on the scales of our interest 
(scales much larger than the Jeans length of baryons). Equation 11411 can be simplihed by integrating it by parts: 
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A further simplification can be made by using an approximation to the integral of the product of spherical Bessel 
functions for Z ^ 1: 

2 7°° 
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Fig. 1.— Simplified schematic diagram illustrating the nature of the correlation between the Doppler and 21-cm anisotropies. Red arrows 
pointing away from the observer indicate ionized gas falling into the positive density perturbation (represented by the black oval) from 
the near side, whereas blue arrows represent ionized gas falling in from the far side. During reionization, there is more ionized gas on the 
near side of the perturbation (at lower redshift) than on the far side. This implies that the net effect from this perturbation is a redshift 
of the CMB in that direction (labeled as < 0)- Because the sources responsible for reionization are located in halos which are very 

biased relative to the underlying linear density field, the overdense region shown here is actually underdense in neutral hydrogen, so that 
this overdensity represents a negative fluctuation in the 21-cm signal (labeled as ^ 21 —cm < 0). Because both the 21-cm and the Doppler 
fluctuations from a region that is undergoing reionization are both the same sign, the signature of reionization is a positive correlation, 
while recombination (in which the situation is reversed for the Doppler signal) results in an anti-correlation. In reality, the growth of 
fluctuations and the dependence of the density on redshift complicate the picture, so that the sign of the signal is determined not by the 
derivative of the ionized fraction d[xe\/dz, but rather d[xe{z){l z)^^'^]ldz (see equation 1 ^ . 


where r(z) = rjo — ri(z) is the comoving distance out to an object at a given z. We obtain 








re 


(17) 


In what follows we will use the exact expression given by equation 03 in our main quantitative results, while we will 
retain the approximate expression given by equation 113 to develop a more intuitive understanding of the origin of 
the cross-correlation. We have found that the exact expression gives results which are about 10% lower (at I ~ 100) 
than the approximate expression of equation 03 for the single reionization history we will use in § 3.4, while for 
the double reionization history the exact result is smaller by about 40%. This is because the line-of-sight integral 
in equation da acts to smooth out features in redshift, an effect which dissappears when the delta function is used 
in the approximation. Since the double reionization model fluctuates much more strongly in redshift than the single 
reionization model, the effect is more apparent for double reionization. 

Equation 03 implies one important fact: the cross-correlation vanishes if Dre is constant. In other words, the 
amplitude of the signal directly depends on how rapidly structure grows and reionization proceeds, and the sign of 
the correlation depends on the direction of reionization (whether the universe recombines or reionizes). Moreover, the 
shape of PCi{z) directly traces the shape of the matter power spectrum at fc = Ijriz'). It is well known that P(fc) has a 
broad peak at the scale of the horizon size at the epoch of matter-radiation equality, fceq — 0.011 Mpc“^ (flm^^/O.lS). 
Since the conformal distance (which is the same as the comoving angular diameter distance in flat geometry) is on the 
order of 10^ Mpe at high redshifts, the correlation power spectrum will have a peak at degree scales, I ~ 10^. 


3.2. Ionized Fraction-Density Correlation 

While Pss{k) is a known function on the scales of interest here, the cross-correlation between ionized fraction and 
density, Pxs{k), is not. In order to understand its importance in determining the observable signal, we have estimated 
its value on large scales. Since we give the full details of derivations in Appendix, we quote only the result here: 

Xe{z)Pxs{k) = -xh{z) \axH{z) \J>h{z) - 1 - /] Pss{k), (18) 


where bh{z) is the average bias of dark matter halos more massive than m„ 




bh{z) = 1 -f 


TT fcol\{z)D{z) <^min 


(19) 
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Tn-mm is the minimum halo mass capable of hosting ionizing sources, /coii(-z) is the fraction of matter in the universe 
collapsed into halos with m > mmin, and (Tmin = cr(mmin) is the r.m.s. of density fluctuations at the scale of mmin 
at 0 = 0. We take mmin to be the mass of a halo with a virial temperature Tmin, which we will treat as a free 
parameter. Here, / is a parameter characterizing the physics of reionization: / = 0 is the “photon counting limit”, in 
which recombinations are not important in determining the extent of ionized regions. On the other hand, / = 1 is the 
“Stromgren limit”, in which the ionization rate is balanced by the recombination rate, as would occur in a Stromgren 
sphere. While our choice for the range of / is resonable, larger values are possible if recombinations limit the size of 
H II regions and the clumping factor increases with increasing density. Equation m is general in the sense that it can 
accomodate such a scenario. It is easy to check that P^s naturally satisfies the physical constraints: it vanishes when 
the universe is either fully neutral, xr = 1, or fully ionized, xr = 0. Although we have derived an explicit relationship 
between PxS and Pss (which is based on several simplifying assumptions - see Appendix), we note that the formulae 
presented here are sufficiently flexible so that any model for the large-scale bias of reionization with PxS = bxsPss can 
be substituted for the one we use here. 

We simplify equations Cl and d by making the approximations that e ” « 1 (justified by observations of the 
CMB polarization; see Kogut et al. 2003) and 


D{z) 


1 -|- Zr 
1-bz ’ 


( 20 ) 


which is a very good approximation at z ^ 1 when the universe is still matter-dominated. The linear growth factor 
has been normalized such that D(zn) = 1; thus, 


d I + Zr _ (Dmgo)^^^(l + Zn) 
dz 1 + Z (1-1- z)l/2 


( 21 ) 


We also use the relation 


f(z) =CrT —(1 - Tp)(l -I- z)'^Xe{z) 

rup 

= 0.0525Hof^6h(l -f zfxe{z), (22) 

where pw is the baryon density at present, Yp = 0.24 is the helium mass abundance (hydrogen ionization only is 
assumed), and Xe(z) is the ionized fraction. One finds 


— {Die ’’) = — 0.0525i?gOmDf,/i(l-I- zn)(1 + Xe{z){l + 


Combining equations d and 121, we obtain 


(23) 
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/o lO^Mpc 

The approximation of equation d and equation (El imply 


ji [Hvo - v{z))]ji [fc(?7o - v{z'))] ■ 


(24) 
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PCf^-^{z) 

27T 


18.4 pK^ 



X 


Pss[l/r{z), zr]{1 + zr)'^ 
10^ Mpc^ 



{z){l + zf/^ 



(26) 


Note that P{k,ZR){l + zrY is independent of the normalization epoch, zr, for zr ^ 1; the result is independent 
of the choice of zr, as expected. These equations are the main result of this paper, and we shall use these results 
to investigate the properties of the correlation in more detail. Since we have a product of dxe/dz and xr, we expect 
that the largest contribution comes from the “epoch of reionization” when Xe{z) changes most rapidly. Therefore, 
by detecting the Doppler-21-cm correlation peak(s), one can determine the epoch(s) of reionization. The sign of the 
cross-correlation is also very important. The sign of the cross-correlation is determined by the sign of the derivative 
term and the difference between xrPss and XePxS- For the case in which xrPss > XePxS, the Doppler effect and the 
21-cm emission are anti-correlated when the ionized fraction, Xe{z), increases toward low z faster than (1 -|- z)~^P. 
For our simplihed model of inhomogeneous reionization (see Appendix), we find that xrPss < XgPxS, however, and in 
this case we find a positive correlation as the universe is being reionized. This is unique information that cannot be 
obtained by present means. See Fig. ^for a schematic diagram which describes the nature of the cross-correlation. 
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Fig. 2.— {Left) Power spectrum of the cross-correlation between the cosmic microwave background anisotropy and the 21-cm line 
fluctuations, I{27r). We assume 2 = 15 for a reionization history given by equation uni with a reionization redshift of Zr = 15 

and duration of Az = 0.5. Note that its shape follows that of the linear matter power spectrum, Pss{^) with k ~ l/r{z), where r{z) is the 
comoving angular diameter distance. (Right) Evolution of the peak amplitude of /{27t) at I ~ 100 from a homogeneous reionization 

history described by equation US. for Az = 0.5 and different reionization redshifts, Zr = 7,11,15 and 19, from left to right. 


3.3. Illustration: Homogeneous Reionization Limit 

It may be instructive to study the nature of the signal by taking the “homogeneous reionization limit”, in which the 
ionized fraction is uniform, Sx = 0. Such a situation may be more relevant than our model for biased reionization if, 
for example, the photons responsible for reionization have a very long mean free path or clumping in denser regions 
cancels the effect of bias. In the homogeneous limit, the approximate formula (eq. I2bll implies 


l^Cf^-^iz) Pss[l/r{z),ZN]{l + ZNf 

2^ ^24.5mK 


XXh{z)— Xe{z){l + zf^"^ 



(27) 


One may estimate the amplitude of the signal at the epoch of reionization, z = Zr^ using a duration of reionization 
at Zr, Az as follows (omitting factors of order unity): 


^{Zr) 195 XHiZr)Xe{Zr) / 1 + Zj-\ 

^ “ Ai 0)25 10 ) 


(28) 


The remarkable feature is that the predicted signal is rather large. For Zr = 15 (which is consistent with early 
reionization suggested bv iKogut et alJ l|2003ll 'l and Az = 1, we predict /{2n) ^ —600 at I ^ 10^. 

Under these assumptions, therefore, detection of the anti-correlation peak should not be too difficult, given that the 
Wilkinson Microwav e Anisotropy Probe (WMAP) has already obtained an accurate CMB temperature map at I ~ 10^ 
(jRennett et a.1.ll2003ll . When any experiment for measuring the 21-cm background at degree scales becomes on-line, 
one should correlate the 21-cm data on degree scales with the WMAP temperature map to search for this peak. Note 
also that in the homogeneous reionization limit the sign is reversed, so that reionization results in an anti-correlation. 
The sign of the correlation therefore depends sensitively on the degree to which reionization is biased on large scales. 


3.4. Reionization History 

To calculate the actual cross-correlation power spectrum, we need to specify the evolution of the ionized fraction, 
Xe(z). Computing Xe(z) from first principles is admittedly very difficult, and this is one of the most challenging tasks 
in cosmology today. To illustrate how the cross-correlation power spectrum changes for different reionization scenarios, 
therefore, we explore two simple parameterizations of the reionization history. 

In one case, we ass ume that the ionized fraction increases mono tonically toward low z. We use the simple parame¬ 
terization adopted bv IZaldarriaga. Furlanetto. fc Hernouistl ll2004tl : 
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xh{z) 


1 -I- exp [—(z — Zr)/Az] ’ 


(29) 
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Fig. 3. — {top and left panels) Peak correlation amplitude vs. redshift. Each panel is labeled with a different value of / which 
parameterizes the uncertainty in the physics of reionization (see Appendix for the definition of / and detailed discussion). The most likely 
value of / is somewhere between 0 and 1 . T he dotted line corresponds to the homogeneous reionization limit in which fluctuat ions in the 
ionized fraction are totally ignored (Eq. I27l h while the thick line takes into account fluctuations in the ionized fraction (Eq. 1261 1. The 
dashed line is the difference between the homogeneous reionization and the total signal, {bottom right) Evolution of Te with redshift. Note 
that in all cases the reionization of the universe results in a positive correlation. 


where Zr is the “epoch of reionization” when xnizr) = 1/2 and Az corresponds to its duration. In this case, one 
obtains a fully analytic formula for the correlation power spectrum: 


27T 


~58^K^ [4/3 +lna;//(z) (6/, - / - 1)] 


P[l/r{z),ZN]{l + zn^ f ^bh'^\ f0,mh 


XXHiz)Xe{z) 


3 _ xh{z){1 + z) 
2~ Az 
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10® Mpc^ 

3/2 


0.02 


0.15 


1/2 


(30) 


In the homogeneous reionization limit, PxS = 0, one gets PCi/{2tt) ~ —165 for z = 9 = Zr and Az = 1/2, and 
the amplitude of the signal scales as (1 + z^)®/^, as expected (see Eq. EHl). For an early reionization at Zr = 15, the 
homogeneous reionization model predicts PCi/{2tt) ~ —570 
The left panel of Figure [3 shows the absolute value of the predicted correlation power spectrum, Z^C'f^“^/(27r), for 
the homogeneous reionization model with z = 15 = z^ and Az = 0.5. As we have explained previously, the shape 
of P\CP~^\ exactly traces that of the underlying linear matter power spectrum, Pss- The right panel of Figure [3 
shows the the redshift evolution of the peak value of the power spectrum at I ^ 100, for different values of z^. As 
discussed at the end of § 3.2, the reionization of the universe leads to an anti-correlation between the Doppler and 
21-cm fluctuations. The magnitude of the signal increases with redshift when the duration of reionization in redshift, 
Az, is fixed (see equation ®)- We could instead hx the duration of reionization in time. At, in which case Az 
increases with redshift as Az oc (1 + z)®/^At; according to equation (EHI . therefore, the peak height in this case would 
be approximately independent of redshift. 

To gain more insight into how the prediction changes with the details of the reionization process, let us use a 
somewhat more physically motivated model for the ionized fraction, 

ln[l - Xe(z)] = -Coiz)fco\liz). 


(31) 
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Fig. 4.— Same as in Fig. 2, but for a “double reionization” model in which the universe undergoes a brief period of recombination. 
Note that in all cases the recombination epoch results in a negative correlation. 


The ionized fraction increases monotonically toward low 2 when ("o does not depend on 2 . Using this model with 
(0 = 200 and Tmin = lO"* K, we calculate the cross-correlation power spectrum. Figure |21 plots the peak value of 
as a function of z, showing the contribution from Pgs, PxS, and the sum of the two (Eq. I^^). The bottom- 
right panel shows the evolution of the ionized fraction predicted by equation m- In this figure we explore the 
dependence of the signal on the details of reionization by varying the parameter /. (See Appendix for the precise 
meaning of /.) In all cases, the contribution from Pss is negative, whereas that from P^s is positive; because the halo 
bias is relatively large for our fiducial case of Tmin = 10"* K, with A <bh < 17 for 10 < z < 30, the P^s term dominates 
over the Pss term, and the correlation is positive (see also Iliev et al. 2005). Increasing the value of / towards a more 
recombination dominated scenario decreases the importance of the dominant PxS term, reducing the total amplitude 
of the signal further. What happens when the universe was reionized twice (Cen 2003; see however Furlanetto & 
Loeb 2005)? In Figure 01 we showed the case where the ionized fraction is a monotonic function of redshift. As seen 
in the figure, there is a prominent correlation peak, regardless of the details of reionization process, encoded in /. 
The situation changes completely when the universe was reionized twice. We parameterize such a double reionization 
scenario using a z-dependence for rmin(z) and Co(-z): 


and 

where 


Co(2) = Ci + (C/ - Ci)9{z) 

Tminiz) = Ti + {Tf - Ti)g{z), 

\ _ exp [ (Zcrit '2^)/Aztran] 

1 T exp [ (Zcrit z')/ Aztran] 


(32) 

(33) 

(34) 


is a function that approaches zero for z > Zcrit and unity for z < Zcrit) with a transition of duration Aztran- We take 
Zcrit = 15, Aztran = 0.25, = 100, ()/ = 40, Ti = 10^ K, and Tf = 10"^ K. In this case, the minimum source halo 

virial temperature makes a smooth transition from 10^ K at high redshift to 10"* K at low redshift, as might occur if 
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dissociating radiation suppresses star formation in “minihalos” with virial temperatures < 10^ K. The drop in Co('2^)j 
which for convenience coincides with the transition in Tmin, could be due to, for example, metal pollution from Pop 
III stars creating a transition to Pop II, accompanied by a transition from a very top heavy IMF to a less top heavy 
one (e.g. Haiman & Holder 2003). In this scenario, the universe may recombine until enough Pop II stars and halos 
with virial temperatures > lO^iF form to finish reionization. We emphasize that this is a simple parameterization 
for illustration, and is not meant to represent a realistic double reionization model. However, this model is sufficient 
to show that the signature of a recombination epoch during reionization is a reversal in the sign of the correlation. 
Because of a rapid change in the ionized fraction during recombination, the negative correlation peak is very prominent, 
reaching ~ 700 — 900 for / = 0 — 1. 


4. PROSPECTS FOR DETECTION 
4.1. Error Estimation 


Assuming that CMB and instrumental noise for 21-cm lines are Gaussian, one can estimate the error of the correlation 
power spectrum by 


where Al is the size of bins within which the power spectrum data are averaged over I — Al/2<1<1 + Al/2, and /sky 
is a fraction of sky covered by observations. 


/sky = ^ = 2.424 X 10-3 
47r 


n 


100 deg^ 


(36) 


In the I range we are considering {I ~ 10^), CMB is totally dominated by signal (i.e., noise is negligible), which gives 
^ 2 (^cmb/( 27 r) ^ (50 /iK)^ at I ~ 10^. On the other hand, the 21-cm lines will most likely be totally dominated by 
noise and/or foreground and the intrinsic signal contribution to the error may be ignored. We also assume that the 
fo reground cleaning reduces it to below the nois e level. We calculate the noise power spectrum based upon equation (59) 
of llZaldarriaga. Furlanetto. fc HernauistI II2004H : 


27r Autohs V 27r A/T) ’ 


(37) 


where Av is the bandwidth, fobs is the total integration time, and A/T is “sensitivity” (an effective area divided by 
system temperature) measured in units of m^ The maximum multipole, Zmax, for a given baseline length, D, is 

given by 


^max = 2tt^ = 2994 


D 

1 km 


10 

1 + z 


(38) 


Here, we have used A = 21 cm(l -|- z). Note that we have implicitly assumed uniform coverage for the interferometer 
in deriving equation ezj, which may not be realistic. Making the baseline distribution more compact would enhance 
the detectability. 


4.2. Square Kilometer Array 

The current design of the Square Kilometer Array (SKA) aims at a sensitivity of A/T ~ 5000 m^ K-^ at 200 MHz. ^ 
Of which, 20% of the area forms a compact array configuration within a 1 km diameter, whereas 50% is within 5 km, 
and 75% is within 150 km. Since we are interested in a relatively low-l part of the spectrum, we use the compact 
configuration, D = 1 km, and A/T = 5000 x 0.2 = 1000 m^ K-^. We obtain 

_ (130 AtK)2 

27r A%oiith AiqviHz 

where A%onth is the number of months of observations and Az^muz is the bandwidth in units of MHz. Note that we 
have assumed here that all the time spent during the observation is on-source integration time. However, a more 
realistic assesment would be that a smaller (e.g., ^ 1/3) fraction of the total time is spent integrating. In this case, 
one should compensate by increasing the total time of the observation accordingly. Since the noise power spectrum is 
much larger than the amplitude of the predicted correlation signal, we safely ignore the contribution of Cf^~^ to the 
error. (We ignore the second term on the right hand side of Eq. 1351 b 

The planned contiguous imaging field of view of SKA is currently 1 deg^ at A = 21 cm and it scales as A^. 
Using the number of independent survey helds, Ageid, the total solid angle covered by observations is given by D ~ 
100 deg^[(l -I- z)/10]^lVfieid. This estimate is, however, based on the current specification for the high frequency 
observations, and may not be relavant to low frequency observations that we discuss here. It is likely that there will 

® Information on SKA is available at http://www.skatelescope.org/pages/concept.htm 


I 

100 


1 + z 
10 


D 

1 km 


10^ m^ K- 
A/T 


(39) 
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be different telescopes with much larger field of view at low frequencies, and thus we shall adopt a field of view which 
is three times larger: 


~ 300 deg^ 



(40) 


Using equation ll^ and the parameters 
order of 


Err 



of SKA, we find the expected error per VA^monthAfieidA ^tvihz to be on the 


~ 938 


l/Al 


;2(7cmb/(27r) 


A^monthA^field^^MHz 2500 ^K^ 


(41) 


Therefore, for the nominal survey parameters, A^month = 12 and A^fieid = 4, the SKA sensitivity to the cross-correlation 
power spectrum reaches Err[Z^C';^^“^/(27r)] ~ 135 /rK^, which gives ~ 3-(T detection of the correlation peak for 
the normal reionization model, and ~ fi-cr detection of the anti-correlation peak for the double reionization model. 
Increasing the integration time or the survey fields will obviously increase the signal-to-noise ratio as VA^monthA^fieid- 
One would obtain more signal-to-noise by choosing a larger value for Az/, which is equivalent to stacking different 
frequencies. (However, Av must not exceed the width of the signal in frequency space.) Therefore, we conclude that 
the cross-correlation between the Doppler and 21-cm fluctuations is fairly easy with the current SKA design. For 
more accurate measurements of the shape of the spectrum, however, a larger contiguous imaging field of view may be 
required. The more promising way to reduce errors may be to increase sensitivity (i.e., larger AjT) by having more 
area, A, for the compact configuration. This is probably the most economical way to improve the signal-to-noise ratio, 
as the error is linearly proportional to {A/T)~^, rather than the square-root. 


5. DISCUSSION AND CONCLUSION 

We have studied the cross-correlation between the CMB temperature anisotropy and the 21-cm background. The 
cross-correlation occurs via the peculiar velocity field of ionized baryons, which gives the Doppler anisotropy in CMB, 
coupled to density fluctuations of neutral hydrogen, which cause 21-cm line fluctuations. Since we are concerned with 
anisotropies in the cross-correlation on degree angular scales {I ~ 100), which correspond to hundreds of comoving 
Mpc at z ^ 10, we are able to treat density and velocity fluctuations in the linear regime. This greatly simplifies the 
analysis, and disting uishes our work from previous work o n similar subjects that dealt only with the cross-correlation 
on very small scales llCnnra,vll2flfl4t ISalvaterra, et a, 1.1120(1511 . Furthermore, because the 21-cm signal contains redshift 
information, the cross-correlation is not susceptible to the line of sight cancellation that is typically associated with 
the Doppler effect. Finally, because the systematic errors of the 21-cm and CMB observations are uncorrelated, the 
cross-correlation will be immune to many of the pitfalls associated with observing the high redshift universe in 21-cm 
emission, such as contamination by foregrounds^. We argue that detection of the predicted cross-correlation signal 
provides the strongest confirmation that the signal detected in the 21-cm data is of cosmological origin. Without using 
the cross-correlation, it would be quite challenging to convincingly show that the detected signal does not come from 
other contaminating sources. 

We find that the evolution of the cross-correlation with redshift can constrain the history of reionization in a 
distinctive way. In particular, we predict that a universe undergoing reionization results in a positive cross-correlation 
at those redshifts, whereas a recombining universe results in a negative correlation (this dependends on our simplified 
model of biased reionization - a model in which reionization is homogeneous would imply a reversal of the sign of 
the correlation). Thus, the correlation promises to reveal whether the universe underwent a period of recombination 
during the reionization process (e.g., Cen 2003), and to reveal the nature of the sources of ionizing radiation responsible 
for reionization. The signal we predict, on the order of PCi/{2 tt) ~ 500 — 1000 /iK^, should be easily detectable by 
correlating existing CMB maps, such as those produced by the WMAP experiment, with maps produced by upcoming 
observations of the 21-cm background with the Square Kilometer Array (SKA). 

Our derivation of the cross-correlation rests upon linear perturbation theory and the reasonable assumption that 
the sizes of ionized regions are much smaller than scales corresponding to I ~ 100. However, assuming that the sizes 
of ionized regions are much smaller than the fluctuations responsible for the signal we predict is not equivalent to 
assuming that the ionized fraction is uniform. Because our prediction depends on the correlation between ionized 
fraction and density, P^s, we have derived a simple approximate model for it (see Appendix). In future work, we 
will use large-scale simulations of reionization to verify the accuracy of the relation we derive, and perhaps to refine 
our analytical predictions. Whatever the result of more detailed future calculations, we are confident that the CMB 
Doppler-21-cm correlation will open a new window into the high redshift universe and shed light on the end of the 
cosmic dark ages. 


We would like to thank Asantha Cooray and the anonymous referee for useful comments. E. K. acknowledges support 
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^ A potential source of foreground contamination is the Galactic synchrotron emission affecting both the CMB and 21-cm fluctuation 
maps; however, the amplitude of synchrotron emission in the CMB map at degree scales is much smaller than the Doppler anisotropy from 
reionization, and thus it is not likely to be a signiflcant source of contamination. 
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APPENDIX 

DENSITY-IONIZATION CROSS-CORRELATION 

The size of H II regions during reionization is a function of the neutral fraction: as the neutral fraction decreases, 
the typical size of H II regions increases, quickly approaching infinity as the neutral fraction approaches zero and the 
H II regions percolate. As predicted by analytical and numerical studies of the large scale topology of reionization 
(Furianetto, Zaldarriaga, & Hernquist 2005; Iliev et al. 2005), the typical H II region size approaches only up to a few 
tens of comoving Mpc at even near the end of reionization. Since we are interested in epochs during which the ionized 
fraction is about a half, we can safely assume for our purposes that the typical H II region size is smaller than the 
length scales of the fluctuations relevant here 100 Mpc). In this case, the ionized fraction within a given volume 
can be determined by considering only sources located inside that volume. 

Let us suppose that we take a region in the universe which has an overdensity of S, where d <C 1. If we assume that 
each baryon within a collapsed object can ionize C(^) baryons, then the ionized fraction within some volume can be 
written as a function of its overdensity S, 

ln[l - a:e((5)] = -C{S)fco\i{S), (Al) 

where /coii(d) is a local fraction of the collapsed mass to mean mass density, which would be different from the average 
collapsed fraction in the universe, /coii(O). Note that this functional form correctly captures the behavior at low and 
high ionizing photon to atom ratio. For C/coii 1, Xe — C/coii 'C 1, which corresponds to all the ionizing photons 
emitted within the volume ionizing atoms within that volume, as expected before H II regions have percolated. For 
C/coii 2> 1, which corresponds to many more ionizing photons than atoms, Xe — 1, as expected after percolation. Given 
that we are only considering sources located within the region, however, this expression is only an approximation during 
percolation, when sources from outside of the volume become visible. We emphasize that in any case equation EH is 
based on a simplifying assumption and does not capture many of the subtleties included in more sophisticated models 
of reionization. 

Here, we present two functional forms for C((5) which are meant to bracket two important physical limits. The first 
limit we will refer to as the “Stromgren limit”, while the second we will refer to as the “photon counting limit”. In 
both limits, we will assume that each hydrogen atom in a collapsed object will produce £ 7 ( 2 ) ionizing photons. 

Stromgren Limit 

If we assume that a fraction r]^.{z) of collapsed gas is undergoing a burst of star formation of duration At^:{z), then 
the ionizing photon luminosity per unit volume, A5y, is given by 

iV, = (A2) 
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where n is the mean density of the universe. The “Stromgren limit” is defined such that every recombination is 
balanced by an emitted photon; the following equation therefore applies: 


ln[l - Xe{5)\ 


Nr,, 


aciAt*n(l +15) ’ 


(A3) 


where a is the recombination coefficient, c/ is the clumping factor, and Nr,, is the recombination rate per unit volume 
in a fully-ioNaed IGM. The last two terms in equation iXdli are the ratio of photon luminosity within a given volume 
to the number of recombinations per unit time which would occur in that volume were it to be fully ionized. Note 
again that this expression ensures the proper behavior of x, in the low and high photon luminosity limits. Combining 
equations EH and (IA3II . we find 

C(<5) =Co(l + (5)-'~Co(l-<5), (A4) 

where (^o = e 7 ? 7 */(ac;ACri) and the approximation is valid in the limit 5 <?; 1. In deriving this relation, we have 
assumed that the clumping factor, ci , is independent of 5. While it is unlikely that clumping decreases with increasing 
5, it is plausible that it could increase. This would decrease the correlation between density and ionized fraction, 
PxS- Because this term typically dominates the cross-correlation (see § 3.4), this would have the effect of reducing the 
predicted signal. 


Photon Counting Limit 

In the “photon counting limit”, recombinations are not important in determining the the extent of ionized regions. 
Instead, it is the ratio of all ionizing photons ever emitted to hydrogen atoms which determines the ionized fraction. 
The number density of photons that have been emitted within a volume with overdensity 5 is given by 

n~^{5) = £ 771(1 -b (5)/coii((5). (A5) 

In the photon counting limit, we assume that the ionized fraction is given by 

In [I - Xe{S)] = = _g (A6) 

77(1 -b 0 ) 

Combining equations EH and (IA5II . we find that ( is independent of overdensity, 

C(<5) = Co = £ 7 - (A7) 


Dependence of Collapsed Fraction on S 
Motivated by these two limits, we parameterize the ^-dependence of C as 

C(^) = Co(l-/<5). 


(AS) 


In the photon counting limit, / = 0, while in the recombination dominated limit / = I. Now that we have specified 
the form of C(^); we turn to the collapsed fraction, /coii(^). The average collapsed fraction in the universe, i.e., /coii(<5) 
with (5 = 0, is given by 

/coii(O) = erfc ^ 


According to the extended Press-Schechter theory, the local collapsed fraction is (Lacey & Cole 1993) 

Sc{z) - S/D{z) 


fcoiiiS, m) = erfc 




(A9) 


(AIO) 


where m is the mass of the region. On large scales, a{m) CTmin and 5 I, so that equation (rnm can be expanded 
in a Taylor series around i5 = 0, 

fcoiiiS) ~ /coii(O) -b J - . . (5. (All) 

V TT aminD(z) 

An alternative expression for the local collapsed fraction can be written in terms of the mean halo bias bh, 


fconiS) = /con(0)i-^, 

I -b 0 

which, for (5 <C I, is well-approximated by 

fcollid) ~ /coll(O) [l -b (6;t — I)i5] . 

Equations (TniTl and im are consistent only if 




(AI2) 

(AI3) 


bh 1 + 


TT /coll(O) 


(AM) 
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Fig. A5.— The effect of bias on the relative importance of the and terms of the cross-correlation signal. The dotted, solid, and 
dashed curves in each panel correspond to models in which the sources responsible for reionization have a minimum virial temperature of 
10^, 10"^, and 10^ K, respectively, {bottom-right) Shown are the reionization histories with a value of ^0 chosen such that the Thomson 
scattering optical depth Tes = 0 . 15 . The model with Tmin = 10 "^ K is a reionization history obtained with (^0 = 200 , and is the same as 
the single-reionization mode l pres ented in the main body of the paper, (bottom-left) Mean halo bias vs. redshift. (top-left) Ratio Pxs/Pss 
as obtained using equation iA20i . (top-right) The ratio \xePx5]/t^H Pss]’ Note that when the ratio equals one, the correlation vanishes, 
while larger values indicate a dominant contribution from the P^-s term. The top panels assume the photon counting limit (/ = 0 ). 


The reader can easily verify that this expression is the same as that found by averaeing over the halo bias derived by 
Mo & White (1996), 


v'^ — 1 
= 1 + c . 

(A15) 

which gives 

T fZ. 

(A16) 

where 

f{iy) (X exp [-i^V2] . 

(A17) 

The Taylor expansion of the collapsed fraction used in deriving equation lAllll is 
standard linear bias formalism. 

therefore fully consistent with the 

Final Expression 

Equations lAlll. (IA8II. and 1A13II imply that 


ln[l - Xe{5)] = ln[l - Xe] [l + (bh - 1 - f)S] 

(A18) 


where we have used the fact that ln[l — Xe] = — Co/coii(0). In the limit where / = 0 and we are simply counting 
photons, the ionized fraction does not depend on 6 if the mean halo bias = 1, since the additional photons emitted 
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within a given region are exactly canceled by the additional atoms contained within that region. If recombinations are 
important, however, the condition for Xe to remain constant is given by = 2. In this case, the bias must compensate 
for the additional photons necessary to balance the enhanced recombination rate per atom within the volume. As was 
noted in § A.l, this relies on assuming a clumping factor which is independent of density. If the clumping factor is 
an increasing function of density, then the condition would he bh > 2. The cross-correlation of ionized fraction and 
density fluctuations is given by (for (5^1 and (5^; <C 1) 

(5x5) ln(l - Xe)(bh - 1 - /)(55), (A19) 

Xe 

SO that 

XePxsik) ~ -(1 - Xe) ln(l - Xe)(bh “ 1 “ f )Pss{k) 

= -XHlnxH(bh - i - f)Pss{k). (A20) 

A comparison of the different terms implied by equa tion (IA20II is shown in Figure ^5] 

If the universe begins to recombine, then equation would be correct in the limit where the recombination time 

is short compared to the time it tak es for sources to dininish in intensity. For simplicity, we will assume this is the 
case and use the relation of equation exclusively in the main body of the paper. In the following section, we will 

investigate the departure from that relation for the case where the sources decay faster than the recombination time. 


Bias in a Recombining Universe 

Equation was derived under the assumption that the ionized fraction is determined by the abundance of 

reionization sources and the density of their environment. However, in the limit where the intensity of ionizing 
radiation due to these sources drops precipitously, as may be expected from metal enrichment or some other form of 
negative feedback, the ionized fraction will be determined by the rate of recombination. In order to understand the 
effect of a “recombination epoch” on the cross-correlation, we will derive a simple relation for P^s for the extreme case 
in which a region of the universe recombines with no sources present. 

In the absence of ionizing radiation, recombination is expected to proceed according to 

^=-(I + 5)x^, (A2I) 

dy 


where y = tjtrec is time in units of the mean recombination of the universe, tree- We will take the initial ionized fraction 
to be a deterministic function of the overdensity, so that the initial fluctuation of ionized fraction (when recombination 
begins occur) is Sx^i = b^^iS, where the subscrip t i re fers to the initial value. If the bias in ionized fraction just before 
recombination begins is described by equation IjAlflII . then we have 


bx,i = ln(I - Xe,i){bh - I - /). 

Xe,i 


(A22) 


For the sake of generality, however, we will report our results in terms of b^^i- Solving for equation dMU, we obtain 
the time evolution of Xe, 

^eiy) = , , "'T. , ^ (A23) 


I -I- Xe,iil -f 6)y’ 


from which it follows that 


^x(Xe) — 

where we have assumed 5 <C I and have used the relation 


^ (bx,^ + 1 ) - 1 

^e.i 


^iv) = 


1 + Xe,iy 


(A24) 


(A25) 


When Xe = Xe,i, Sx = bx,i5, as expected. As the universe recombines to become fully neutral, Xe ^ h and 5x —5. 
Since we expect bx,i > I because overdense regions have an overabundance of ionizing dources, a period of recombination 
is expected to weaken the importance of PxS term. When Xelxe^i = 1/2 and bx,i S> I, for example, the bias determined 
from is too large by a factor of two. Because we have assumed that sources turn off instantaneously in equation 

ifXTlTl . this is an upper limit to the effect of a recombination epoch on PxS- 








